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Abstract. 

In the framework of the theory of open systems based on completely positive 
quantum dynamical semigroups, we give a description of continuous variable quantum 
entanglement and quantum discord for a system consisting of two non-interacting non- 
resonant bosonic modes embedded in a thermal environment. We study the time 
evolution of logarithmic negativity, which characterizes the degree of entanglement, 
and show that in the case of an entangled initial squeezed thermal state, entanglement 
suppression takes place for all temperatures of the environment, including zero 
temperature. We analyze the time evolution of the Gaussian quantum discord, which 
is a measure of all quantum correlations in the bipartite state, including entanglement, 
and show that discord decays asymptotically in time under the effect of the thermal 
bath. We describe also the time evolution of classical correlations and quantum mutual 
information, which measures the total correlations of the quantum system. 

PACS numbers: 03.65.Yz, 03.67.Bg, 03.67.Mn 
1. Introduction 

In recent years there is an increasing interest in using non-classical entangled states 
of continuous variable systems in applications of quantum information processing, 
communication and computation pQ. In this respect, Gaussian states, in particular two- 
mode Gaussian states, play a key role since they can be easily created and controlled 
experimentally. Due to the unavoidable interaction with the environment, in order to 
describe realistically quantum information processes it is necessary to take decoherence 
and dissipation into consideration. Decoherence and dynamics of quantum entanglement 
in continuous variable open systems have been intensively studied in the last years 

In this paper we study, in the framework of the theory of open systems based on 
completely positive quantum dynamical semigroups, the dynamics of continuous variable 
quantum entanglement and quantum discord of a subsystem consisting of two uncoupled 
bosonic modes (harmonic oscillators) interacting with a common thermal environment. 
We are interested in discussing the correlation effect of the environment, therefore we 
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assume that the two modes are independent, i.e. they do not interact directly. The 
initial state of the subsystem is taken of Gaussian form and the evolution under the 
quantum dynamical semigroup assures the preservation in time of the Gaussian form of 
the state. In Refs. [TH [191 120] we studied the evolution of entanglement and quantum 
discord of two identical harmonic oscillators interacting with a thermal environment for 
initial symmetric squeezed vacuum and squeezed thermal states of the subsystem. In 
this work we extend the previous analysis to the case of non- resonant bosonic modes and 
take non-symmetric squeezed thermal states as initial states. We show that in the case 
of an entangled initial squeezed thermal state, entanglement suppression (entanglement 
sudden death) takes place for all temperatures of the environment, including zero 
temperature. We analyze the time evolution of Gaussian quantum discord, which is 
a measure of all quantum correlations in the bipartite state, including entanglement, 
and show that discord decays asymptotically in time under the effect of the thermal 
bath. Before the suppression of the entanglement, the qualitative evolution of discord is 
very similar to that of the entanglement. We describe also the time evolution of classical 
correlations and quantum mutual information, which measures the total correlations of 
quantum system. 



2. Equations of motion for two modes interacting with an environment 



We study the dynamics of a subsystem composed of two non-interacting bosonic modes 
in weak interaction with a thermal environment. In the axiomatic formalism based on 
completely positive quantum dynamical semigroups, the Markovian irreversible time 
evolution of an open system is described by the Kossakowski-Lindblad master equation 
[2"T| 122] • We are interested in the set of Gaussian states, therefore we introduce such 
quantum dynamical semigroups that preserve this set during time evolution of the 
system. The Hamiltonian of the two uncoupled non-resonant harmonic oscillators of 
identical mass m and frequencies uj\ and U2 is 

H=^{pl+p 2 y ) + ^W + ujly\ (1) 

where x, y are the coordinates and p x , p y are the momenta of the two quantum oscillators. 

The equations of motion for the quantum correlations of the canonical observables 
x, y and p x ,p y are the following (T denotes the transposed matrix) [22J: 

da(t) 



dt 



Ya(t) + a(t)Y l +2D, 



(2) 
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quantities. We introduced the following 4x4 bimodal covariance matrix: 
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where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual reduced one-mode states, while the matrix C 
contains the cross-correlations between modes. The elements of the covariance matrix 
are defined as a,^ =< RiRj + RjRi > /2, i,j = 1, .., 4, with R = {x,p x , y,p y }, which up 
to local displacements fully characterize any Gaussian state of a bipartite system. 
The time-dependent solution of Eq. ^ is given by [22] 

cr(t) = M(t)[a(0) - <x(oo)]M T (t) + a(oo), (5) 

where the matrix M(t) = exp(Yt) has to fulfill the condition lim^oo M(t) = 0. The 
values at infinity are obtained from the equation 

Ya(oo) + a(oo)Y T = -2D. (6) 



3. Dynamics of quantum correlations 

3.1. Time evolution of entanglement 

In order to quantify the degree of entanglement of the two-mode states it is appropriate 
to use the logarithmic negativity. For a Gaussian density operator, the logarithmic 
negativity is completely defined by the symplectic spectrum of the partial transpose of 
the covariance matrix. It is given by En = max{0, — log 2 2z>_}, where z/_ is the smallest 
of the two symplectic eigenvalues of the partial transpose a of the two-mode covariance 
matrix a [6]: 



2v- 



r,2 



A =f V A 2 - 4deta 



, - ■ v - (7) 

and A is the symplectic invariant (seralian), given by A = det A + det B — 2 det C. 
In our model, the logarithmic negativity is calculated as 



E N {t) = max{0, ~ \og 2 [4g(a(t))]}, 



(8) 



where 




(det A + det B) - det C 



(det A + det B) - det C 



1/2 



det a(t) 



(9) 



It determines the strength of entanglement for Eif{t) > 0, and if £?jv(i) = 0, then the 
state is separable. 
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We assume that the initial Gaussian state is a two-mode squeezed thermal state, 
with the covariance matrix of the form [251 
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(10) 



with the matrix elements given by 

a = n\ cosh 2 r + n<i sinh 2 r + - cosh2r, (11) 

b = ni sinh 2 r + n 2 cosh 2 r + - cosh 2r, (12) 

c = ^(^i +^2 + 1) sinh2r, (13) 

where ni,n 2 are the average number of thermal photons associated with the two modes 
and r denotes the squeezing parameter. In the particular case n% — and n 2 = 0, 
(10) becomes the covariance matrix of the two-mode squeezed vacuum state. A two- 
mode squeezed thermal state is entangled when the squeezing parameter r satisfies the 
inequality r > r s [23], where 

cosh2rs = K±iK^±i) (14) 
m + n 2 + i v ; 

We suppose that the asymptotic state of the considered open system is a Gibbs state 
corresponding to two independent bosonic modes in thermal equilibrium at temperature 
T. Then the quantum diffusion coefficients have the following form (we put from now 
on h = 1) [21]: 



muiD xx = PxPx = - coth 



mui 2 2kT 

D X p x — D yPy — D xy — D PxPy — D xpy — D yPx — 0. 

The evolution of entangled initial squeezed thermal states with the covariance 
matrix given by Eq. (10) is illustrated in Fig. 1, where we represent the dependence 
of the logarithmic negativity -EW(t) on time t and temperature T for the case of an 
initial non-symmetric Gaussian state (a ^ b). For all temperatures T, including zero 
temperature, at certain finite moment of time, which depends on T, E^{t) becomes 
zero and therefore the state becomes separable. This is the so-called phenomenon of 
entanglement sudden death. It is in contrast to the quantum decoherence, during which 
the loss of quantum coherence is usually gradual [T2J One can also show that 

the dissipation favors the phenomenon of entanglement sudden death - with increasing 
the dissipation parameter A, the entanglement suppression happens earlier. The same 
qualitative behaviour of the time evolution of entanglement was obtained previously 
[T9t [27] in the particular case ri\ = and n 2 = corresponding to an initial two-mode 
squeezed vacuum state and in the case of symmetric initial squeezed thermal states. 
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Figure 1. Logarithmic negativity En versus time t and temperature T for an 
entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 3, rii — 3, — 1 and A = 0.1, toi = 1, ui2 = 2. We take m = % = k = 1. 

Comparing the present results with those obtained in the previous paper [20] 
one can assert that the asymmetry (a ^ b) of the initial Gaussian state favors the 
suppression of entanglement. The most robust under the influence of the environment 
is the entanglement of symmetric (a = b) initial squeezed thermal states. An even 
stronger influence on the entanglement has the non-resonant character of the two modes: 
by increasing the ratio of the frequencies of the two modes, the entanglement sudden 
death happens earlier in time. The longest surviving entanglement takes place when 
the modes are resonant (oj\ =002). This effect due to the non- resonance of the modes 
is stronger for small values of the frequencies, and it diminishes, for the same ratio of 
frequencies, by increasing the values of frequencies. 

In our model, in which we suppose that the asymptotic state of the considered open 
system is a Gibbs state corresponding to two independent bosonic modes in thermal 
equilibrium, an separable initial state remains separable in time, and it is not possible 
to generate entanglement. This is in contrast with the possibility of entanglement 
generation starting, for instance, with a separable state in the case of two non-interacting 
two- level systems immersed in a common bath [28]. At the same time we remind that 
we have studied previously [2"ll |2"9] the evolution of the entanglement of two identical 
harmonic oscillators interacting with a general environment, characterized by general 
diffusion and dissipation coefficients, and we obtained that, for separable initial states 
and for definite values of these coefficients, entanglement generation or a periodic 
generation and collapse of entanglement take place. In discussing the entanglement 
decay, it is interesting to mention that models have been elaborated to realize quantum 
feedback control of continuous variable entanglement for a system consisting of two 
interacting bosonic modes plunged into an environment, based on a local technique [30J, 
or on a nonlocal homodyne measurement [31] . 
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3.2. Gaussian quantum discord 

Recent studies have shown that separable states, usually considered as being classically 
correlated, might also contain quantum correlations. Quantum discord was introduced 
[32] 133] as a measure of all quantum correlations in a bipartite state, including - but 
not restricted to - entanglement. Quantum discord has been defined as the difference 
between two quantum analogues of classically equivalent expressions of the mutual 
information, which is a measure of total correlations in a quantum state. For an arbitrary 
bipartite state p\ 2 , the total correlations are expressed by quantum mutual information 

/(Pi 2 )= E S( Pi )-S(p 12 ), (16) 

i=l,2 

where pi represents the reduced density matrix of subsystem i and S(p) = — Tr(plnp) 
is the von Neumann entropy. A measure of bipartite classical correlations C(p\ 2 ) in the 
bipartite quantum state p\2 based on a complete set of local projectors {IT^} on the 
subsystem 2 can be given by 

C(pi 2 ) = S(Pi) - inf {n fe } {5(pi| 2 )}, (17) 

where S(pi\ 2 ) = Y^kP k S(.Pi) is the conditional entropy of subsystem 1 and inf {S(pi\ 2 )} 
represents the minimal value of the entropy with respect to a complete set of local 
measurements {nJ;}. Here, p k is the measurement probability for the kth local projector 
and p\ denotes the reduced state of subsystem 1 after the local measurements. Then 
the quantum discord is defined by 

D(p 12 )=I(p 12 )-C(p 12 ). (18) 

Originally the quantum discord was defined and evaluated mainly for finite 
dimensional systems. Recently [MJ ES] the notion of discord has been extended to 
the domain of continuous variable systems, in particular to the analysis of bipartite 
systems described by two-mode Gaussian states. Closed formulas have been derived for 
bipartite thermal squeezed states [31] and for all two-mode Gaussian states [35J. 

The Gaussian quantum discord of a general two-mode Gaussian state p 12 can be 
defined as the quantum discord where the conditional entropy is restricted to generalized 
Gaussian positive operator valued measurements (POVM) on the mode 2, and in terms 
of symplectic invariants it is given by (the symmetry between the two modes 1 and 2 is 
broken) [35] 

D = /(ft) - f{P-) - f(v+) + /(Vi), (19) 



where 



_ , , OC I 1 _ OC I 1 1 _ iXj 1 , _ _ « 

/(*) = —z- log — n 7T- log —z-, (20) 
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Figure 2. Gaussian quantum discord D versus time t and temperature T for an 
entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 3, Tlx = 3, ri2 = 1 and A = 0.1, uj\ = 1, u)2 = 2. We take m = K = k = 1. 



£ = < 
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(21) 
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a 



4 det A, /3 = 4det5, 7 = 4detC, 5 = 16 det a, 
and u T are the symplectic eigenvalues of the state, given by 



2i£ 



A =f VA 2 - 4 det a, 



:p _ , v _ (23) 

where A = det A + det B + 2 det C. 

The evolution of the Gaussian quantum discord D is illustrated in Figure 2, where 
we represent the dependence of D on time t and temperature T for an entangled initial 
non-symmetric Gaussian state, taken of the form of a two-mode squeezed thermal state 



(10), for such values of the parameters which satisfy for all times the first condition in 



formula (21). The Gaussian discord has nonzero values for all finite times and this fact 



certifies the existence of non-classical correlations in two-mode Gaussian states, either 
separable or entangled. Gaussian discord asymptotically decreases in time, compared 
to the case of logarithmic negativity, which has an evolution leading to a sudden 
suppression of entanglement. For entangled initial states the Gaussian discord remains 
strictly positive in time and in the limit of infinite time it tends asymptotically to zero, 
corresponding to the thermal product (separable) state, with no correlations at all. 

From Fig. 2 we notice that, in agreement with the general properties of the 
Gaussian quantum discord [35], the states can be either separable or entangled for 
D < 1 and all the states above the threshold D = 1 are entangled. We also notice 
that the decay of quantum discord is stronger when the temperature T is increasing. 
It should be remarked that the decay of quantum discord is very similar to that of the 
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Figure 3. Degree of classical correlations C versus time t and temperature T for 
an entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 3, n\ — 3, rii — 1 and A = 0.1, uji — 1, 102 — 2. We take m = h = k = 1. 



entanglement before the time of the sudden death of entanglement. Near the threshold 
of zero logarithmic negativity (E^ = 0), the nonzero values of the discord can quantify 
the non-classical correlations for separable mixed states and one considers that this fact 
could make possible some tasks in quantum computation [36J. The discord is increasing 
with the squeezing parameter r and it is decreasing with increasing the ratio of the 
frequencies u\ and U2 of the two modes and the difference of parameters a and b. 



3.3. Classical corellations and quantum mutual information 

The measure of classical correlations for a general two-mode Gaussian state P12 can also 
be calculated and it is given by [35] 

c = /(V^)-/(V£), (24) 

while the expression of the quantum mutual information, which measures the total 
correlations, is given by 

/ = f{Va) + /(#) - f{v-) ~ f{"+). (25) 

In Figs. 3 and 4 we illustrate the evolution of classical correlations C and, 
respectively, quantum mutual information J, as functions of time t and temperature 
T for an entangled initial Gaussian state, taken of the form of a two-mode squeezed 



thermal state (10). These two quantities manifest a qualitative behaviour similar to 
that one of the Gaussian discord: they have nonzero values for all finite times and in 
the limit of infinite time they tend asymptotically to zero, corresponding to the thermal 
product (separable) state, with no correlations at all. One can also see that the classical 
correlations and quantum mutual information decrease with increasing the temperature 
of the thermal bath. One can show that the classical correlations and quantum mutual 
information increase with increasing the squeezing parameter r and the difference of 
parameters a and b. At the same time classical correlations increase with the ratio 
of the frequencies U\ and u 2 of the two modes, while quantum mutual information is 
decreasing with increasing this ratio. 
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Figure 4. Quantum mutual information / versus time t and temperature T for 
an entangled initial non-symmetric squeezed thermal state with squeezing parameter 
r = 3, m = 3, ri2 = 1 and A = 0.1, =l,u)2 = 2. We take m = h = k = 1. There are 
also represented the Gaussian quantum discord and classical correlations. 



For comparison all these quantities are represented also on the same graphic in 
Fig. 4. In the considered case the value of classical correlations is larger than that of 
quantum correlations, represented by the Gaussian quantum discord. 

4. Summary 

We investigated the Markovian dynamics of quantum correlations for a subsystem 
composed of two non-interacting bosonic modes embedded in a thermal bath. We have 
analyzed the influence of the environment on the dynamics of quantum entanglement 
and quantum discord for Gaussian initial states. We have described the time evolution of 
the logarithmic negativity in terms of the covariance matrix for non-symmetric squeezed 
thermal states, for the case when the asymptotic state of the considered open system is a 
Gibbs state corresponding to two independent quantum harmonic oscillators in thermal 
equilibrium. The dynamics of the quantum entanglement strongly depends on the 
initial states and the parameters characterizing the environment (dissipation coefficient 
and temperature). For an entangled initial squeezed thermal state, entanglement 
suppression (entanglement sudden death) takes place for all values of the temperatures 
of the environment, including zero temperature. The time when the entanglement is 
suppressed decreases with increasing the temperature and dissipation. 

We described also the time evolution of Gaussian quantum discord, which is a 
measure of all quantum correlations in the bipartite state, including entanglement. 
The values of quantum discord decrease asymptotically in time. This is in contrast 
to the sudden death of entanglement. The time evolution of quantum discord is very 
similar to that of entanglement before the sudden suppression of the entanglement. 
Quantum discord is decreasing with increasing the temperature. After the sudden 
death of entanglement the nonzero values of discord manifest the existence of quantum 
correlations for separable mixed states. We described also the time evolution of classical 
correlations and quantum mutual information, which measures the total correlations of 
the quantum system. 
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